Statistics of Bubble Rearrangements in a Slowly Sheared Two-dimensional Foam 
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Many physical systems exhibit plastic flow when subjected to slow steady shear. A unified picture 
of plastic flow is still lacking; however, there is an emerging theoretical understanding of such flows 
based on irreversible motions of the constituent "particles" of the material. Depending on the 
specific system, various irreversible events have been studied, such as Tl events in foam and shear 
transformation zones (STZ's) in amorphous solids. This paper presents an experimental study of the 
Tl events in a model, two-dimensional foam: bubble rafts. In particular, I report on the connection 
between the distribution of Tl events and the behavior of the average stress and average velocity 
profiles during both the initial elastic response of the bubble raft and the subsequent plastic flow at 
sufficiently high strains. 

PACS numbers: 82.70.-y,83.60.La,62.20.Fe 



I. INTRODUCTION 

Bubble rafts have been used as a model experimen- 
tal system for the study of crystalline and amorphous 
solids P, and for the study of two-dimensional foam 
0, 01 . This overlap is just one of many examples that 
points to an important question in the study of the me- 
chanical response of materials. Under conditions of slow 
steady shear, what, if any, is the connection between the 
response of "mesoscopic" materials, such as foams, emul- 
sions, pastes, and slurries, and plastic flow of "molecular" 
systems, such as amorphous solids? Based on macro- 
scopic measurements, the systems are similar. There is 
an initial elastic response for small strains (or stresses) 
and a yield stress, above which irreversible deformations, 
or plastic deformations, occur. Eventually, above some 
critical stress (or strain), the system enters a "flowing" 
state that is characterized by irregular periods of stress 
increase and decrease. This is often referred to as un- 
bounded plastic flow. For the purposes of this paper, 
this will simply be referred to as plastic flow. As one 
reduces the shear rate, the critical stress approaches the 
yield stress in such a way that for sufficiently slow shear 
rates the behavior of the system is essentially shear-rate 
independent. This is often referred to as the quasistatic 
regime. A complete "microscopic" picture of plastic flow 
still does not exist, where microscopic refers to the fun- 
damental length scale relevant to the system in question. 
For example, in bubble rafts, it would be the dynamics of 
individual bubbles. Open questions include: the micro- 
scopic source of the stress release events; the spatial and 
temporal distribution of such events; and the nature of 
such events during periods of stress increase. Experimen- 
tally, the challenge is identifying systems for which the 
microscopic events are directly observable. This is one of 
the main advantages of mesoscopic systems, such as the 
bubble raft, and the reason for the interest in making 
connections between mesoscopic systems and molecular 
systems, such as amorphous solids. 

Models and simulations of diverse systems, ran ging 
from solids 0, 0,111111113,111111 to foam 



IHIIIIIIIIIIIIEIEIIII, have provided a number 
of insights into these questions. The focus of this paper 
is on the role of Tl events in foam. Aqueous foam con- 
sists of gas bubbles separated by liquid walls jH 0, |2l| . 
A Tl event is an irreversible neighbor switching event. 
Both a schematic representation and an actual Tl event 
are presented in Fig. For the purposes of this paper, 
we will focus on the role of Tl events during the steady 
shear of foam. However, it should be mentioned that 
understanding the nonlinear events that are not shear 
induced may be important when comparing foam and 
amorphous solids. In the absence of external stress, foam 
coarsens, and Tl events occur due to geometric changes 
in the foam structure. These Tl events are not necessar- 
ily distinguishable from those caused by flow. In contrast, 
most amorphous solids do not exhibit the equivalent of 
coarsening. However, thermally activated events may be 
important. Possible differences between thermally acti- 
vated and coarsening events is an interesting question, 
but one beyond the scope of this paper. 

Simulations of flowing foam have characterized differ- 
ent aspects of Tl events. Often, one separately consid- 
ers Tl events and "avalanches", i.e. sudden releases of 
stress (or energy) in the foam. One issue is whether 
or not the probability of the number of Tl events in a 
given avalanche exhibits power-law behavior. A common 
model to study dry foam (where the bubbles are essen- 
tially polygonal) is the vertex model J^] . In the 
vertex model, a Tl event is defined to occur when the 
distance between two vertices (i.e. the wall between two 
bubbles) is below a threshold value. One then elimi- 
nates this wall, creating a Tl event. Therefore, within 
this model, all Tl events are essentially instantaneous. 
In this case, simulations found evidence for power-law 
behavior of the probability of Tl events. Another char- 
acterization of the Tl events is the number of Tl events 
per bubble per unit strain, Rti- For the vertex model, 
this quantity is Rti — 0.5. A modified version of the 
vertex model has been used to study the issue of flow 
localization under shear 1^^- These simulations report a 
correlation between the spatial localization of Tl events 
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to the neighborhood of a system boundary with the lo- 
cahzation of shear in the same region. The issue of shear 
locahzation will be discussed in more detail later. 

Another class of models focuses on wet foam (foam 
in which the bubbles are essentially spherical, or in 
two-dimensions, circular), using a quasistatic simulation 
These simulations involve making a small step 
strain and then allowing the system to relax to an energy 
minimum before applying the next step strain. Anytime 
the energy decreases after a step strain, one declares this 
an "avalanche" or "event" . By comparing neighbors in 
the initial and final state, one can count the number of Tl 
events for a particular avalanche. In these simulations, 
avalanches consisting of a large number of Tl events were 
observed, suggesting the possibility of power-law behav- 
ior For this model, Rti was not reported. 

Wet foam under steady shear has also been simulated 
using a q-Potts model |23|. In this case, different bubbles 
are identified by different spin orientations. Simulations 
of the q-Potts model find that the distribution of topo- 
logical rearrangements are not power-law-like. However, 
the distribution of energy drops may be consistent with 
power-law behavior |2fl| | . 

Another important set of simulations for wet foam in- 
volved studying the steady shear of the bubble model 
[T^ IrR l2ll |. This model treats bubbles as spheres (or 
circles) that interact via a spring force proportional to 
their overlap and a viscous drag proportional to velocity 
differences. As this model directly simulates dynamics, 
the duration of Tl events exhibits a distribution of du- 
ration times. Simulations focused on small shear-rates in 
the quasistatic limit, i.e. the flow properties were in- 
dependent of the shear rate. Under these conditions, 
no evidence of power-law behavior is observed in the 
bubble model at high bubble density, and Rti = 0.15 
[TgL HtL l2ll |. If one decreases the density of the bubbles, 
it appears that the distribution of events approaches a 
power-law as one approaches the critical density for the 
"melting" of the foam 21j. 

Before discussing the current state of experiments in 
foam, it is useful to put the theoretical work on Tl events 
in foam in the context of two points of view of plasticity 
in amorphous materials. First, the idea of shear transfor- 
mation zones (STZ's), as developed by Falk and Langer 
0, has received significant attention. STZ are a way 
of describing local, irreversible rearrangements of parti- 
cles during shear. STZ are based on previous work by 
Spaepen and Argon on activated transitions and Turn- 
ball, Cohen, and others on free-volume fluctuations. As 
the STZ refers to a small region of the material with cer- 
tain properties |26j |. there is only a loose connection be- 
tween the STZ and Tl events. However, it is reasonable 
to identify as an STZ regions in which a few Tl events 
combine to form a local slip (see, for instance. Fig. O;). 
It is expected that the local rearrangements identified 
as STZ are associated with quadrupolar energy fluctua- 
tions. In fact, the expected quadrupolar energy fluctua- 
tions have been observed associated with Tl events in a 



simulation of foam |22|, but not, as of yet, in simulations 
of molecular systems. 

Another view of plasticity is based on shear induced 
changes in the potential energy landscape, as proposed 
by Malandro and Lacks 0|. This picture derives from 
an inherent structure formalism and focuses on changes 
in the macroscopic mechanical response of a material due 
to shear induced changes of the potential energy. This 
formulism has been used to study simulations of a qua- 
sistatic version of the bubble model [2^ . In this case, sys- 
tem wide rearrangement events are observed. This is not 
seen in bubble model simulations of the quasistatic limit, 
but it is seen in other quasistatic simulations of foam. 
The work in Ref. suggests the need to carefully de- 
fine the concept of an "event" , especially for steady-state 
experiments where the time scale for events to occur rela- 
tive to the applied shear can be important. For example, 
a shear-rate regime may exist that is quasistatic as de- 
fined by the behavior of quantities such as the average 
stress, but not in a quasistatic with regard to the dura- 
tion of stress releases. Hence, large events get "broken 
up" by the steady shear, changing the nature of the dis- 
tribution of events. 

A number of experimental studies of bubble rearrange- 
ments in model foam have been carried out. As men- 
tioned, some of the earliest work was done using bubble 
rafts 0,0123 1 i-S- layers of gas bubbles floating on a liq- 
uid surface, as a model molecular system, both for crys- 
talline and amorphous solids One major advan- 
tage of bubble rafts is that their two dimensional nature 
allows for easy imaging and tracking of all of the "par- 
ticles" in the system. Another reason that bubble rafts 
have been so useful in the study of molecular systems is 
that there exists detailed calculations of the bubble in- 
teractions [S^l- More recently, bubble rafts were used to 
study rearrangements after a step strain in order to make 
comparisons with the quasistatic simulations of foam |^ . 
This work did not directly measure Tl events, but it did 
look at changes in the number of neighbors for bubbles. 
The results suggested that large scale events were possi- 
ble. 

Experiments have also been carried out using mono- 
layer foam |3l| . Langmuir monolayers consist of a sin- 
gle layer of molecules confined to the air-water inter- 
face. They exhibit a large number of two-dimensional 
phases, including gas-liquid coexistence. This allows for 
the formation of a foam of gas bubbles with liquid walls. 
For a monolayer foam under steady shear, only a small 
number of simultaneous Tl events were observed, with 
Rti — 0.12 ± 0.03. These results are consistent with the 
bubble model. 

As mentioned, the other aspect of Tl dynamics is their 
potential role in explaining shear localization in yield- 
stress materials 22], such as foam and granular systems. 
It has long been known that a yield stress and/or nonlin- 
ear viscosity can lead to inhomogeneous flows . How- 
ever, it is only recently that experimental techniques have 
allowed for detailed measurements of such behavior. A 
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number of such studies have been carried out in granular 
materials, where exponential velocity profiles (or other 
strongly localized velocity profiles) are generally observed 
|33l I34L Issf . In contrast, measurements in various three- 
dimensional pastes, slurries, and foams show a different 
type of inhomogeneous flow. In this case, the flow is not 
strongly localized, and there is a shear-discontinuity at 
the boundary between flow and no flow |36l Is^ . 

For two-dimensional foams, the situation is ambiguous. 
Three-dimensional foam that is confined between plates 
to form a model two-dimensional system exhibits shear 
localization analogous to granular systems This 
work motivated simulations of the modified vertex model 
discussed earlier that showed a connection between TI 
events and shear localization 22]. In this case, it ap- 
pears that the spatial distribution of stress released by 
the TI event results in the subsequent localization of the 
events. The localization of TI events is correlated with 
the localization of flow. In contrast, experiments with a 
bubble raft exhibit a shear-discontinuity 1391 similiar to 
that reported in Refs. js^ls^l- In Ref. [sj, TI events 
were not measured. 

The work reported in this paper addresses the gen- 
eral question of the temporal and spatial distribution of 
TI events during the slow, steady shear of a bubble raft. 
Also, connections between the TI events and the velocity 
profiles reported in Ref. [s^l are made. The rest of the 
paper is organized as follows. Section ^provides the de- 
tails of the experimental setup. The results are presented 
in two parts. Section IlII Al presents the initial response 
of the system. Section Fill Bl presents the behavior dur- 
ing plastic flow. Finally, the results are summarized and 
discussed in Sec. IIVI 



II. EXPERIMENTAL METHODS 

The experimental system consisted of a standard bub- 
ble raft 01 in a Couette geometry. The bubble raft was 
produced by flowing regulated nitrogen gas through a hy- 
podermic needle into a homogeneous solution of 80% by 
volume deionized water, 15% by volume glycerine, and 
5.0% by volume Miracle Bubbles (Imperial Toy Corp.). 
The bubble size was dependent on the nitrogen flow rate, 
which was varied using a needle valve. A random distri- 
bution of bubble sizes was used, with an average radius 
of I mm. The resulting bubbles were spooned into a 
cylindrical Couette viscometer. This produced a two- 
dimensional model of a wet foam on a homogeneous liq- 
uid substrate. Figure |5] presents a schematic side view of 
the bubbles in the apparatus and an image of a top view 
of the bubble raft. 

Due to the nature of the bubble raft, no measurable 
coarsening was observed. However, after approximately 
two hours significant numbers of bubbles would pop, pre- 
sumably due to loss of fluid in the walls. This set the 
upper limit on the total time of the measurements. In 
contrast, during the initial two hour period, only two out 



of approximately 400 bubbles in the field of view were 
observed to pop. 

An important feature of the bubble raft is the gas area 
fraction. To achieve a desired gas area fraction, the bub- 
ble raft was constructed by placing the approximate num- 
ber of desired bubbles in the trough with the outer bar- 
rier set to a large radius. It is important to note that the 
bubbles exhibited a strong attraction to each other. The 
outer barrier was compressed until the desired radius was 
reached. The gas area fraction was determined by thresh- 
olding images of the bubbles and counting the area inside 
of the bubbles. Because of the three-dimensional nature 
of the bubbles, this represents an operational definition 
of gas-area fraction based on the details of the image 
analysis. However, the choice of threshold was consistent 
with an estimate of the gas area fraction based on the 
area of trough and expected distribution of bubble sizes. 
For all of the data reported here, the gas area fraction 
was approximately 0.95. 

The Couette viscometer is described in detail in 
Ref. li^l and shown schematically in Fig. [S^a). It con- 
sists of a shallow dish that contains the liquid substrate. 
Two concentric Teflon barriers are placed vertically in 
the dish. Sections of both of these barriers are visible in 
Fig. EJb) . The outer barrier is a ring consisting of twelve 
segmented pieces. It has an adjustable radius. For the 
experiments discussed here, the outer radius was fixed at 
To — 7 A3 cm. The inner barrier, or rotor, is a Teflon 
disk with a radius = 3.84 cm. The outer edge of the 
disk is a knife edge that is just in contact with the water 
surface. It was suspended by a wire to form a torsion 
pendulum. 

To shear the bubble raft, the outer Teflon barrier was 
rotated at a constant angular velocity Q = 8x 10"'' rad/s. 
The first layer of bubbles at either boundary did not slip 
relative to the boundary. Due to the finite size of the bub- 
bles, this results in an effective inner radius of r = 4.4 cm. 
Due to the cylindrical geometry, the shear rate is not uni- 
form across the system and is given by 7(r) = J' 7^-7^- 
Here v{r) is the azimuthal velocity of the bubbles. During 
plastic flow, the average azimuthal velocity of the inner 
cylinder is zero. Measurements of the average azimuthal 
velocity profile allows for calculations of the shear rate. 
As measured at r = 4.4 cm, 7 = 4 x 10~^ s~^. In this 
regime, where reported, the strain (7) is taken to be the 
strain at this radius, and is computed from 7 = 7i -f 70, 
where t is the time since the initiation of plastic flow and 
7o is the amount of strain developed during the initial 
period. During the initial period, the inner barrier has 
a finite angular speed. However, one can still compute 
the effective shear rate at r = 4.4 cm. In this regime, 
7 = 3 X 10~^ s~^. Again, where reported, the strain is 
the strain at the inner cylinder: 7 = jt, where in this 
case t is measured from the initiation of shear. 

The details of the velocity measurements are given in 
Ref [33 • Video tape of roughly one third of the trough 
was recorded. Images from this tape were taken every 
3.2 s and digitized. An image processing routine based 
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on standard Labwindows functions was developed that 
detected and tracked individual bubbles. This tracking 
software was also used to compute the average bubble dis- 
placements. This is used to compute the deviation of the 
bubble motion from ideal elastic behavior, as discussed 
in Sec. ITTm 

The Tl events were measured by stepping the digi- 
tized images one frame at a time and visually searching 
for the location and time at which Tl events occurred. 
A Tl event was defined to occur when two bubbles were 
observed to lose contact, and two other bubbles moved 
into the resulting space (see Fig.^. Due to the associ- 
ated motions of the other neighboring bubbles, Tl events 
are relatively easy to detect by hand l4l| . For automatic 
tracking of Tl events, it is critical to accurately detect 
essentially all of the bubbles, as one is interested in de- 
termining neighbor switching events. This is in contrast 
with the displacement and velocity measurements where 
the requirement is that one tracks enough bubbles to have 
sufficient statistics. These are the reasons that automatic 
methods were used for displacement and velocity mea- 
surements, but the detection of Tl events was by hand. 

The stress on the inner rotor was determined using two 
different methods. In both cases, the torque, T = kO, 
on the inner rotor is determined by measuring the an- 
gular displacement, 9. (For the experiments presented 
here, the torsion constant k = 5.7 x 10~^ Nm.) The 
stress is then determined from a = T/27rr^. The differ- 
ence in the two methods is the determination of 9. The 
first method uses a magnetic flux technique, and the de- 
tails of this technique are in Ref. This is the more 
precise of the two methods, with a stress resolution of 
3 X 10~'^ mN/m. The second method uses the video im- 
ages of the inner cylinder and tracks fixed features on the 
disk. This method has a resolution of 0.043 mN/m. The 
second method is used to correlate the video analysis of 
bubble motions (displacements and Tl events) with the 
detailed stress fluctuations determined from the magnetic 
flux measurements that are reported in Ref. 42] . 

As mentioned, foams are inherently nonequilibrium 
systems. One complication that arises from this is the 
definition of the yield stress. For sufficiently low shear 
rates, foam will spontaneously release stress, usually as 
part of coarsening process as bubble sizes change. This 
complication is minimized in the bubble raft given that 
no substantial coarsening was observed in the absence of 
applied shear. In either useful operational def- 

inition of the yield stress is the zero shear-rate limit of 
the stress. For the bubble raft of interest here, the aver- 
age stress as a function of shear rate is well described 
by a Herschel-Bulkley model ia{j) ^ To + nj") Q. 
From these results, one can determine a yield stress: 
To = 0.8 ± 0.1 mN/m jH El- For the particular shear 
rate of interest here, this is different from the "critical" 
stress at which the system begins to undergo "steady" 
flow. 



III. EXPERIMENTAL RESULTS 

A. Elastic Regime 

The initial stress response of the system is given in 
Fig. 121 There are a number of interesting features of 
this regime. First, there are three separate regions of the 
initial response, which is essentially set by the slope of 
the stress versus strain curve. These regions are indicated 
by the vertical dashed lines and are separated by isolated 
stress-drops. 

The first region is the linear, elastic response of the 
material. During this period no Tl events are observed. 
The second two regions represent plastic deformations 
in the sense that Tl events occur. These events are too 
small to produce stress drops. But, they modify the slope 
of the stress-strain curve and generate irreversible defor- 
mation. Hence, the identification of these regions with 
plastic deformations. The onset of plastic response is an- 
other useful definition of the yield stress. However, there 
is always ambiguity associated with the definition of the 
onset of Tl events due to the possibility of Tl events 
that are the result of coarsening and not shear. As dis- 
cussed for the bubble raft, this difficulty is minimized 
as no coarsening was observed. However, for the mea- 
surements in Fig. 13 only a fraction of the trough was 
in view. This limits the degree to which the yield stress 
can be measured by this method. However, it is useful to 
note that the onset of Tl events for the single set of data 
studied here is consistent with the value of yield stress as 
determined by fits to the behavior of stress as a function 
of rate of strain. 

The fact that any deviation of the stress-strain curve 
from linear behavior is small allows for the definition of 
an effective shear modulus of the bubble raft, G, during 
periods of stress increase. For the initial region in Fig.O 
G is the elastic shear modulus. The calculation of G as- 
sumes that the stress is proportional to the strain. The 
boundary conditions consist of a fixed rotation rate at 
the outer boundary and an inner boundary that is free 
to rotate, but supported by a torsion wire. Because of the 
symmetry of the Couette geometry, the azimuthal veloc- 
ity, v{r) is only a function of the radial position r. This is 
a standard problem; however, given the slightly unortho- 
dox boundary conditions of this experiment, the solution 
is repeated here. The relevant constitutive equation is 

air) = G7(r). (1) 

Here j{r) is the shear strain, and <j{r) is the resulting 

shear stress. In the cylindrical geometry, 7(r) = ''-^f-^, 
where 9{r) is the angular displacement of the bubble raft. 
For a material confined between two cylinders, the shear 
stress is given by a{r) — r/(27rr^). This follows directly 
from balancing torques on each material element. If the 
bubble raft was a perfectly rigid solid, one would simply 
have v{r) = fir. This is due to the fact that the inner 
boundary is supported by a torsion wire and rotates as 
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it measures the torque on the inner cyUnder. However, 
for a finite value of G, plugging into Eq. ^ we get 



T 



G 



dO{r) 
dr 



(2) 



Integrating this equation, and using the fact that the 
bubble raft does not slip at either boundary, gives for G, 



G = 



and for 



v(r) = fi?' - 



KLUr 






AttG 
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This equation for v can be rewritten by plugging in for 
G, 



vr) = Qr + 



R^ -r 



r:' 



(5) 



The second piece in the expressions for v{r) is due to the 
elastic nature of the bubble raft and the motion of the 
inner cylinder. When n/G is small, the system behaves 
as a rigid body {uj = Q), as expected (large G limit). 

Using the above results, one can find G from measure- 
ments of the average velocity using Eq.0]or from u using 
Eq. 13 For example, the results for v{r) in region (A) of 
Fig. O are given in Fig.0] The solid line is a fit to Eq. 0] 
with v{r)/r = 8.2 x 10^"* rad/s - 0.003 rad/cm2s(l/r2 - 
0.0193 cm~^). The Tl events result in a reduction in 
the effective elastic modulus of the bubble raft. The 
calculated values of G for the three regions are: (A) 
G = 11.2 ± 0.1 mN/m; (B) G = 5.4 ± 0.1 mN/m; and 
(C) G = 8.9 ±0.1 mN/m. 

We use the elastic regime to provide a characterization 
of the local deviation from clastic flow. First, we take the 
flt of v{r) in region (A) as the dcflnition of "ideal" elastic 
motion. Knowing this velocity curve, we can compute the 
expected displacement of a bubble during a strain inter- 
val. From this, we define A to be a measure of the devia- 
tion from elastic behavior: A = ^y{x — Xe)"^ + iu — TJeY, 
where x and y are the actual displacements of the bub- 
ble and Xe and j/e are the expected displacements if the 
motion was "ideal" elastic behavior. As can be seen from 
Fig. 0] even in the "pure elastic" regime there is a sig- 
nificant non-zero variation to the bubble motions. (The 
error bars represent the standard error based on the stan- 
dard deviation of measured velocities in each radial bin.) 
The variation in bubble velocity is due to a combination 
of effects, including the obvious fact that one expects a 
distribution of displacements due to the finite size of the 
bubbles. Therefore, the variation in displacements from 
the ideal elastic behavior in region (A) is used to set a 
minimum threshold value for A of 0.05 cm. Bubbles with 
a value of A below this threshold are considered to have 
undergone "elastic" motion. Even with this cutoff, there 



are a small number of bubbles in the tails of the distribu- 
tion that are classified as deviating from elastic behavior 
even in region A. This is illustrated in Fig.|Sli-c. Each of 
Fig. EJi-c, represent a period of strain of 0.064 in which 
no Tl events occur. The periods were selected from the 
corresponding region (A - C) of Fig. 13 The circles indi- 
cate the location of tracked bubbles (so only a fraction of 
the total bubbles are shown). The color of the bubbles 
indicates the deviation from elastic behavior, with white 
bubbles having a value of A < 0.05 cm. The color code 
is indicated in the figure. 

Figure|3li and e illustrate two classes of Tl events that 
do not result in a stress drop. Fig. is from region (B) 
of Fig. 121 This illustrates an isolated pair of Tl events 
that have an associated region in which the bubbles de- 
viate from elastic behavior. Fig. |3l3 is from region (C) of 
Fig. 13 This illustrates the slippage of two, short rows 
of bubbles due to simultaneous Tl events. Again, there 
is a relatively localized region of deviation from elastic 
behavior associated with these Tl events. 



B. Plastic Flow Regime 

In the plastic flow regime, there are two main ques- 
tions regarding the Tl events. First, what is the corre- 
lation between Tl events and stress? Second, what is 
the correlation between Tl events and bubble motion, as 
measured by either the average velocity or the deviation 
from elastic behavior? 

Regarding the correlation between Tl events and 
stress, it is interesting to consider the periods of stress in- 
crease. As with region B and C in Fig.|2| there are often 
Tl events during these periods of stress increase. There- 
fore, in general, these are periods of plastic deformation, 
though preliminary observations suggest that occasional 
increases exist during which no Tl events occur. One 
way to characterize stress increases is to use Eq.|3to cal- 
culate an effective shear modulus, G, for each separate 
period of stress increase. This can than be correlated 
with the number of Tl events that occur in that period. 
A preliminary measurement of this is shown in Fig. 
This result is preliminary because only a fraction of the 
sample was viewed. Therefore, the results for the num- 
ber of Tl events represent a lower bound. However, it is 
interesting that the data all falls below the straight line, 
suggesting a correlation between G and the number of 
Tl events, as expected. 

It is natural to expect that one necessary condition for 
a Tl event to occur is that the local stress exceeds some 
critical value. This would suggest a correlation between 
the stress and the location of the Tl events. Figure [7| il- 
lustrates that no correlation exists between the stress on 
the inner cylinder and the radial positions of Tl events. 
One would expect such a correlation if the stress field was 
given by the continuum limit, which in a Couette geome- 
try is a{r) = {a{ri)rf)/r'^, and the critical stress for a Tl 
event was spatially uniform. Under these conditions, for 
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each cr(ri), there is a maximum r at which Tl events can 
occur. This is set by the critical stress required for the 
generation of a Tl event. The lack of a correlation sug- 
gests that at least one, if not both, of these assumptions 
is false. In fact, work in other systems suggests that both 
assumptions are false. Given the direct measurement of 
stress chains in granular matter 'SS^, it is reasonable to 
expect such chains in the bubble raft. This would rep- 
resent a breakdown of the continuum assumption for the 
stress distribution. Also, simulations of amorphous metal 
have shown the existence of localized, high stress regions 
(referred to as r-defects) that are the sources of local 
flow |12|. In other contexts, models and simulations have 
su ggest ed the existence of "weak" zones in complex fluids 
p4l lifiL liq that are the source of viscous-like behavior. 
A more detailed study of these issues will required im- 
proved images, but the current work is very suggestive. 

To summarize the average properties of Tl events as a 
function of strain. Fig. |H1 plots the number of Tl events 
per bubble versus strain. Again, this is shown simultane- 
ously with the stress versus strain curve to illustrate the 
general correlation between the size of the stress drops 
and the total number of Tl events. One observes that 
most stress drops consist of a cascade of events through- 
out the stress drop. As with the stress increases, de- 
tailed correlations between the size of a stress drop and 
the number of Tl events will require images of the entire 
sample. However, one can compute Rti- For this shear 
rate, Rti = 0.18 ± 0.01, in reasonable agreement with 
both the bubble model and Langmuir monolayer foam. 

The next question is the connection between velocity 
profiles and Tl events. Based on the results of Ref |39l |. 
it is known that there exists a shear discontinuity at 
Tc — 6.7 cm for the system reported on here. There- 
fore, there is no expectation of strong localization of the 
Tl events are reported in Ref. 22\ because there is no 
shear localization. However, one might expect a connec- 
tion between the radial distribution of Tl events and the 
shear discontinuity. 

The shear discontinuity divides the system into two 
regimes. Below Tc, the average velocity is consistent with 
that of a power-law fluid. Above rc, the systems acts 
like an elastic solid. Figure |51 illustrates the connection 
between the average velocity proflle and the spatial dis- 
tribution of Tl events. The vertical line indicates the 
spatial location of the shear discontinuity [s^ . The basic 
shape of the distribution is similar to that found in the 
simulations reported in Ref. There is a "peak" at 

smaller radii, with the distribution tailing off as one goes 
to larger radii. The main difference is the location of the 
cutoff in the Tl distribution. As reported in Ref. [2^ . 
the cutoff in velocity and Tl events is at essentially the 
same radius. In contrast, for the system reported on 
here, there is no obvious signature in the distribution of 
Tl events at the shear discontinuity (see Fig. This 
may be due to the fact that even though the shear-rate is 
zero, the bubbles are still moving near the outer wall, and 
differences in bubble size may lead to Tl events. Also, 



it may be an artifact of how close the shear discontinu- 
ity is to the outer wall. Future work in larger systems is 
needed to better understand this issue. 

In order to better understand the detailed connection 
between Tl events and stress drops, two short periods 
of strain arc highlighted, as indicated in Fig. ^| These 
are segments of the data presented in Fig. [7| The period 
of strain illustrate in Fig. llOb was selected to highlight 
the nature of potential correlations between the stress 
behavior and the number of Tl events. First, the initial 
elastic rise A is included for comparison with the stress 
increase in the interval labelled C. During the initial rise, 
there is only one observed Tl event, and the effective 
elastic modulus is G = 5.6 ± 0.1 mN/m. In contrast, 
during region C, there are 9 observed Tl events, and the 
effective elastic modulus is G = 2.1 ± 0.1 mN/m. These 
results reinforce the general connection between number 
of Tl events and effective elastic modulus discussed with 
respect to Fig. |S1 In contrast, the regions labelled B - 
E all have roughly the same number of Tl events. Yet, 
region B and E are stress drops. Region C is a period 
of stress increase, and region D is a slight decrease in 
stress. One difficulty in drawing definitive conclusions 
from this data is the fact that only part of the system 
is being viewed. However, this strongly suggests that 
the additional bubble motions, not just the Tl events, 
play an important role in determining the overall stress 
evolution. 

The interval illustrated by Fig. llOb was selected to 
make connections with the velocity profiles reported in 
Ref. [S^ in an attempt to better understand the shear 
discontinuity that occurs at r^.. Here, the average bubble 
displacements are measured, but these are directly re- 
lated to average velocities. This sequence is particularly 
interesting because there are three stress drops that oc- 
cur at different average stress values, and the drop at the 
lowest average stress (region E) exhibits the larger value 
of Tf. [s^. (We are considering the behavior in region B 
and C as two separate stress drops because of the short 
plateau between them. However, this points out the issue 
regarding the definition of "events" as discussed earlier 
in the context of Ref. p^.l 

The spatial distribution of Tl events and bubble de- 
viations from elastic behavior are given in Fig. 1111 using 
the same criteria as described for Fig. O in Sec. 1111 Bl 
White bubbles represent essentially elastic behavior, and 
the color of the other bubbles is the degree to which their 
motion deviates from elastic. The letters in Fig.llOb cor- 
respond to the labelling of Fig. ^2 One observes very 
similar distributions of Tl events for all three stress drops 
(see Fig. [TTb.c. and e). If one looks carefully, the dis- 
tinguishing factor appears to be the number of bubbles 
deviating from elastic behavior at any given radius. This 
is made clearer by considering the average bubble dis- 
placement as a function of radial position, as illustrated 
in Fig. m 

Figure [T^ is a plot of /S.9/Vl/S.t versus radial position. 
The values of are computed by dividing the system 
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into 20 equally spaced radial bins and averaging the an- 
gular displacements over all bubbles in a given bin over 
the time interval of interest. The time intervals are se- 
lected so that they match the strain intervals indicated in 
Fig. llOb . For comparison, the displacements during the 
essentially flat regions in stress are given as open symbols, 
and the displacements during the stress drops are given 
as closed symbols. The angular displacement is normal- 
ized by UAt. The solid line is the "ideal" elastic behavior 
given by the fit to the data in Fig. ^ One can see that 
for interval E (solid triangles), the deviation from the ex- 
pected elastic displacement occurs at the largest value of 
Tc. This is consistent with the results reported in Ref. 
for the velocity profiles. What is new in these results is 
the ability to correlate the location of Tl events during a 
stress drop and the location of the deviation from elastic 
behavior. 

For event E, there are two clear deviations from elastic 
behavior, as shown in Fig. E| (a) at r = 6.48 cm there 
is a positive deviation; and (b) at r = 6.16 cm there is a 
negative deviation. During the drop, the maximal radial 
position of a Tl event is r = 6.24 cm. The existence 
of positive and negative deviations is consistent with the 
bubbles associated with a Tl event moving both forward 
and backward relative to the average flow. The average 
displacements during B and C are essentially identical. 
However, for C more then B there is some evidence for 
a positive and negative deviation at r — 5.84 cm and at 
r = 5.44 cm, respectively. For these drops, the maximal 
radial position of a Tl event is r = 6.29 cm. Comparing 
these numbers strongly suggests that the location of Tl 
events is not the main contribution to the determination 
of the deviation from elastic behavior, and hence, the de- 
termination of re- Instead, it is the detailed motion of 
the surrounding bubbles. Interestingly, the greatest dif- 
ference between the two events in terms of Tl position is 
in the precursor to the drops; yet the precursors have very 
similar angular displacements (open symbols in Fig. I12|l. 
During the interval labelled D in Fig. llOb . one observes 
Tl events as far out as r = 6.98 cm. In contrast, dur- 
ing the interval labelled A in Fig. llOb . one only observes 
Tl events as far out as r = 5.87 cm. Presumably these 
events play an important role in establishing the local 
stress fields that govern the bubble motions during the 
subsequent stress drop. 



IV. DISCUSSION 

Even though various aspects of the work presented here 
are preliminary in the sense that only a portion of the 
entire raft was studied, a number of questions regarding 
the role of Tl events in the macroscopic response of a 
bubble raft to flow have been addressed. First, the con- 
tribution of Tl events to the effective shear modulus was 
considered. Tl events during periods of stress increase 
effectively lower the shear modulus of the bubble raft. 
Similarly, during stress drops, there is a correlation be- 



tween the size of the drop and the total number of Tl 
events. Future work is required to establish a more de- 
tailed correlation between the number of Tl events and 
the effective shear modulus and size of stress drops. 

Correlations between positions of Tl events, average 
stress, individual bubble displacements, and average an- 
gular displacements of bubbles were considered. A gen- 
eral picture that emerges from these measurements is the 
importance of understanding the local stress field and 
the local geometry of bubbles. For example, an investi- 
gation of individual bubble motions before and during a 
stress drop (see discussion of Figs . ITUl and [TT|l shows that 
the radial distribution of Tl events can not be under- 
stood in terms of a simple continuum model and single 
stress threshold. The Tl events in the strain period im- 
mediately prior to a stress drop play an important role 
in establishing the local stress field and geometric rela- 
tions between bubbles that sets the subsequent motions. 
For example, the two different stress drops highlighted 
in Fig. ^1 exhibited similar distributions of Tl events, 
but the deviations from elastic behavior and the average 
displacements were very different. The main differences 
between the events was in the distribution of precursor 
Tl events, not in the average bubble motions. 

The connection between Tl events and the position 
of the shear discontinuity was also considered. Both in 
terms of the average properties (see Fig.|3l and the short 
time motions (see Fig. I12|) . There is no clear evidence 
for a connection between the positions of Tl events and 
the shear discontinuity, but larger system sizes need to 
be studied. However, there may be an indirect connec- 
tion through the stress relaxation and subsequent motion 
of surrounding bubbles. Indirectly, these measurements 
have some potential implications for the simulations of 
the modified vertex model '2^. These simulations illus- 
trate that a localization of Tl events can lead to a shear 
localization. This system does not exhibit localization of 
either the Tl events or the shear. This indirectly sup- 
ports the connection between Tl event localization and 
shear localization. What remains an important questions 
is why would Tl events localize in one case and not the 
other? An obvious difference between the Tl events in 
the bubble raft and in the simulation is the duration of 
the Tl events. In the model, the Tl events all occur on 
a very short time scale, by construction. For the bub- 
ble raft, there is a distribution of times for the duration 
of Tl events. Some events occur very slowly (over 10 - 
20 seconds). It is these difference in duration that may 
modify the impact on the local stress. Again, this point 
to the importance of understanding the local stress fields 
generated by the Tl events, and not just the distribu- 
tion of the events themselves. Furthermore, as part of 
the future work that focuses on local stress fields, it will 
be important to correlate the changes in local stress with 
the duration of the Tl events. 

The measurements reported here focused on bubble 
displacements and Tl events. Where possible, compar- 
isons with the bubble model show quantitative agree- 
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ment, such as for Rti- This adds support to previous re- 
sults with the bubble raft that were also in general agree- 
ment with the bubble model Having strong agree- 
ment between the experiments and a theoretical model 
is useful for the next stage of experimental studies. Es- 
sentially all of the results point to the need for measure- 
ments of the local stress field. Future experimental work 
is planned that will use the bubble distortion as a direct 
measure of local stress, as has been done with other foam 
systems [43, H^, 112, • Close contact with simulations 
that focus on the stress released by Tl events and STZ's, 
as well as experimental studies of granular material, will 
be important for understanding this future work. 
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FIG. 1: (a) Schematic representation of the three main steps 
in a Tl event. The bubbles labelled 1 and 2 are initially 

neighbors. As the bubbles are sheared, all four bubbles meet 
at a vertex. After the event, the bubbles labelled 3 and 4 are 
neighbors, (b) - (d) are three images illustrating an actual 
Tl event in the bubble raft. The location of one Tl event is 
highlighted by artificially coloring the bubbles involved white. 
The images are taken 3.2 s apart and the white scale bar in 
(b) is 2 mm long. 
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FIG. 2: (a) A schematic drawing of the apparatus showing a 
side view. The main elements of the apparatus arc the knife 
edge disk that is supported by a torsion wire and that serves as 
the inner cylinder for the bubbles. There is a separate fixed 
inner cylinder in the fluid (in gray). There is a segmented 
outer cylinder for generating flow, and there is a fixed dish 
that holds the fluid. The bubbles sit on top of the fluid, as 
indicated by the circles, (b) An image from the top of the 
bubbles in the apparatus that shows a portion of both the 
outer and inner cylinder. The black scale bar in the lower left 
corner is 5 mm. 
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FIG. 5: Five images representing typical bubble deviations 
from elastic flow during the initial stress rise. Images (a) - 
(c) are taken from the corresponding three strain intervals 
in Fig. 1^ and show typical strain intervals in which no Tl 
events occur. Image (d) shows a typical localized Tl event 
from region (B) in Fig.|3] Image (e) shows an event composed 
of multiple Tl events in which two rows of bubbles slip pass 
each other. This event is taken from region (C) in Fig. |3] 
The circles indicate the location of a subset of bubbles that 
have been tracked. The sizes of all the circles are the same, 
independent of actual bubble size, for clarity. They are color 
coded based on the deviation from elastic displacements, as 
defined in the text. White represents deviations less then 
0.05 cm. The color bar gives the scale for deviations greater 
than 0.05 cm. The squares represent the location of Tl events, 
where the color equals the time relative to the start of the 
interval. The scale bar in image (a) is 0.5 mm. 



FIG. 4: A plot of v(r)/r versus radial position for the bubble 
motion during interval A in Fig. |3 The solid squares are data 
averaged over all bubbles at a given radial position. The solid 
line is a fit to Eq. 0] 
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FIG. 6: Scatter plot of the number of Tl events during a 
period of stress increase versus the effective elastic modulus 
G for the same period of strain. The solid line is a guide to 
the eye. 




strain 

FIG. 8: The solid line is the same stress versus strain curve 
as shown in Fig. |7| The bars summarize the data in Fig. |7| 
by plotting only the number of Tl events/bubble in a strain 
interval of 0.013. 
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FIG. 7: The individual points are the radial position of Tl 
events as a function of strain. The solid curve is the stress 
as a function of strain (measured using images of the inner 
cylinder) . 



radial postion (cm) 

FIG. 9: The squares are the average velocity of the bubbles 
normalized by Sir as a function of radial position. The bars 
give the total number of Tl events/bubble as a function of 
radial position. The solid lines are guides to the eye. The 
horizontal line is v{r)/Q,r — 1, which corresponds to the mo- 
tion of a rigid body. The vertical line is the location of the 
shear discontinuity, as reported in Ref. |39|l . 
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strain 

FIG. 11: Five images representing typical bubble deviations 
FIG. 10: This is two strain intervals from Fig. 13 showing from elastic flow during the corresponding strain intervals as 
both the location of Tl events (solid squares) and the stress indicated in Fig. EHIb). The circles indicate the location of 
(sohd line, as measured by the magnetic method) as a function a subset of bubbles that have been tracked. The size of the 
of strain. Each interval is further divided into smaller strain circle is the same for aU bubbles for clarity. The circles are 
intervals by the vertical hues. The labels in (b) correspond to color coded based on the deviation from elastic displacements, 
the images in Fig. [TTI as defined in the text. White represents deviations less then 

0.05 cm. The color bar gives the scale for deviations greater 
than 0.05 cm. The squares represent the location of Tl events, 
where the color equals the time relative to the start of the 
interval. The scale bar in image (a) is 0.5 mm. 
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FIG. 12: The average angular displacement of the bubbles 
normalized by QAt as a function of radial position. The dif- 
ferent symbols are for the different strain intervals in Fig. llOb : 
(A) open circles; (B) solid squares; (C) solid circles; (D) open 
squares; and (E) solid triangles. Here At is the time inter- 
val for each strain interval. The solid line is the fit to elastic 
behavior from Fig. |1| 



